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Abstract 

Pm \ Consider a homogeneous time-continuous branching process where 

Mh . individuals have constant birth rate 5, and hfe length distribution Q 

r^ I having mean E{Q) — 1. Let X{u) denote the number of individuals 

jrt ' alive at time u, and assume that X(0) = 1. Let i^ be a positive 

integer and define Ak '■= /g ^{x(u)=K}du, the accumulated time that 
the branching process has exactly K individuals alive. In this paper we 
prove that E{Ak) = 5^~^ / (fc(l V (5)^), irrespective of the life length 
- distribution Q, subject to the normalizing condition E{Q) = 1. 
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^ : 1 Introduction and main result 

en 

Consider a homogeneous time-continuous branching process X having life- 
length distribution Q, where we measure time in the unit of expected life- 
lengths, implying that E[Q) = 1. During their lives, individuals give birth 
^ \ according to independent Poisson processes with constant rate 5, where each 

. . .' birth consists of a single child. Let X{u) denote the number of individuals 

alive at time u, and assume that the process has one ancestor, i.e. that 
X(0) = L 

Such branching processes have been studied extensively, see, for example, 
Jagers [8] and are sometimes referred to as (binary) splitting trees, see, for 
example, Champagnat et al. |3J . It is well-known that the branching process 
is subcritical, critical or supercritical depending on whether 5 is smaller than, 
equal to, or larger than 1, and that P[X{u) — )• oo) > if and only if 5 > 1. 

Let N denote the number of individuals ever born in the branching 
process and let Qi,Q2-, ■ ■ ■ denote the life-lengths of the individuals (for 



example, labelled according to time of birth) . It then holds that 

_ roo ^ /'oo 

y^Qj= x{u)du=y] / ^hxiu) 

~, Jo t^, Jo 



i)=K}du, 



3=1 •'" K=l- 



where all three expressions equal infinity if the branching process grows 
beyond all limits. It is well-known that the expectation of ([1]) is finite if 
and only if 5 < 1, and, using Wald's lemma, that the mean then equals 
E{N)E{Q) = E{N) = l/{l - 6). 

In the current paper we are interested in the expectation of each of the 
terms on the right hand side of ([TJ. The factor K is just a constant, so we 
instead choose to state our result for 
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Ak '■= / ^{x(u)=K}du, and its mean (2) 

Jo 

J hxiu)=K}du] = J P{X{u) = K)du, (3) 

the expected time the branching process has K living individuals, or the 
expected " iT-occupation" time. We have the following result: 

Theorem 1.1 Consider a branching process with birth rate 6 and life-length 
distribution Q, with E{Q) = 1 but otherwise arbitrary. Let X(u) denote the 
number of individuals alive at time u, and assume X{0) = 1. Then, for 
K = 1,2, . . . we have 



/•oo 

E{Ak) = / P(X( 
Jo 



i;K-l 



The proof of Theorem 11.11 is given in Section [2j 

Remark 1. The statement of the theorem was conjectured by Neal [10] 
who also proved it for the case that Q ~ Exp{l) and 6 < 1. In Neal [lOj . 
it was shown that E{Ax) is the key quantity for computing the asymp- 
totic endemic equilibrium distribution for a subcritical branching process 
with immigration and a supercritical, homogeneously mixing SIS epidemic 
model. 

Remark 2. Theorem 1 1.1 1 is an example of an insensitivity result, in that, 
E{Ak) only depends upon Q through its mean. Similar results have been 
observed for queueing networks, see Zachary [13], for an overview. Further- 
more, Theorem 11.11 follows straightforwardly from Zachary [13], Theorem 1 
in the subcritical case 5 < 1. However, Zachary |13| . Theorem 1 does not 
appear to be easy to adapt to critical and supercritical branching processes. 



The total time until extinction is given by T = X^fcLi ^k = /g ^{x{u)>i}du. 
It is well-known that this time has infinite mean in the critical and super- 
critical case {S > 1). An immediate consequence of the theorem tells us 
what the mean equals in the subcritical case. 

Corollary 1.1 Under the same setting as Theorem \l.l\ and for the case that 
6 < 1 we have 

E(T) = f E(A,-) = -'°^<;-^' . (5) 

K=l 

Remark. The case 5 = can be taken as a limit of ([5]) and the mean is hence 
E{T) = 1 as it should. 

Theorem 11.11 and Corollary 11.11 have direct implications for inference on 
6 when data consists of observing Ak or T = Y1k=i ^k (cf- Farrington and 
Grant [5]). 

Corollary 1.2 Let K > 2 and suppose that Ak = tx is observed for the 
branching process described above. If the process is known to be subcritical 
(5 < 1), then the moment estimator of 6 is given by 6 = {KtKY^~^^ ■ If 
the process is super critical the moment estimator is given by 6 = l/{KtK)- 
If instead T = t is observed and the process is known to be subcritical, then 
the moment estimator 6 is the largest solution to 1 — 6 = e~ . 

Remark 1. It might seem unnatural to assume sub- or super-criticality to 
be known when making inference. It is however well-known that branching 
processes possess many non-standard inference features (cf. Guttorp [Gj). 

Remark 2. In terms of applications, branching processes are common 
models for populations but also for epidemics, at least when the outbreaks 
are small. Corollary 11.21 treating estimation problems is hence of interest 
in these situations, see, for example, Farrington and Grant j5j, treating 
inference problems for a related discrete time version of epidemics. 

Remark 3. The transcendental equation 1 — 5 = e~ appears also in 
mathematical epidemiology, where t is a measure of infectiousness of the 
disease and 5 is the ultimate fraction getting infected in the epidemic. From 
this theory (e.g. Diekmann et al. [1]) it is known that there is a unique 
positive solution 5 in (0,1) if and only if t > 1. If t < 1 then the largest 
solution is 5 = 0. 

Knowing that the expected iC-occupation time, E{Ak), is independent 
of the life-length distribution Q, subject to E{Q) = 1, raises the question if 
the result holds even under more general assumptions. The answer to this 



question is negative, at least if we generalize the current model, having con- 
stant birth rate during Q, to a model where the birth rate is inhomogeneous 
(time- varying) together with a random duration of the life-length distribu- 
tion such that the expected number of births still equals 6. For this extended 
model the corresponding i^-occupation time is no longer constant (and the 
same as in Theorem II. ip . A simple illustration of this is for example Q = 1 
and 6 large, e.g. 10, and K = 1 and we compare two time-varying birth 
rates. The first model is where nearly all of the birth intensity comes early 
in life and the second where it comes close to the end of the life period (i.e. 
close to 1). Since the process will probably never return to 1 (being "very" 
super critical), most of the expected time spent with K = 1 individual alive 
comes from before the first birth, which is clearly longer in the second model. 
Also the result no longer holds if we allow the possibility of multiple births 
at each point of Poisson point process with rate 5. It is straightforward to 
construct a counter example comparing Q = 1 and Q ~ Exp(l) since for 
Q = 1 all individuals born at the same time will die at the same time. 

2 Proof of Theorem 11.11 

2.1 Introduction 

The approach we take to prove Theorem 11.11 is to consider phase-type dis- 
tributions (Asmussen et al. [1] ) for Q. Specifically, we take Q to be a finite 
mixture of hypoexponential distributions. That is, we assume that there 
exists m G N such that Q is a mixture distribution oi Qi,Q2, ■ ■ ■ , Qm with 
P{Q = Qi) = Pi i^^iVi = 1) and, for each z, there exists nj G N such that 

Qi ~ Exp(7j,i) + Exp(7i,2) + • • • + Exp(7j,„J (6) 

with YlT=iPiYl^jLilij = 1- Therefore an individual has a lifetime dis- 
tributed according to Qi with probability pi, where the lifetime Qi consists 
of Hi stages each of which lasts an exponential length of time. Thus if 
we know the distribution from which each of the individual lifetimes come 
and the stage at which each individual is at in their lifetime, the branching 
process is Markovian. Any Coxian distribution (Asmussen et al. pj) can 
be expressed as a mixture of hypoexponential distributions, and therefore 
the above class of mixtures of hypoexponential distributions is dense, As- 
mussen et al. [1] . Hence Theorem 11.11 follows trivially if we can show that 
E{Ak) = 5^~^ /{K{1 A 5)^^ holds for the mixtures of hypoexponential 
distributions. 



For i = 1, 2, . . . , m, let kij denote the total number of individuals with 
a lifetime distributed according to Qi who are currently in stage j of their 
lifetime with kj = {ki^i, ki^2, ■ ■ ■ , h,nj and k = (ki, k2, . . . , k^)- To consider 
the transitions to and from state k it is helpful to define Bij to be a vector 
of length X^i^i n-j whose {i,jY^ (X]I=i "-^ + J*'^) element is 1 and all other 
elements are 0. Transitions from state k occur as follows:- 

k - ejj + ejj+i rate hjjij (i = 1, 2, . . . , m; j = 1, 2, . . . , n^ - 1) 

k - Bi^n, rate ki^n,7i,n, (i = 1, 2, . . . , m) (7) 

k + ej 1 late K5pi (i = 1, 2, . . . ,m,), 

with K = '^^lYyiLi^ij- The Markov branching (birth-death) process 
defined by ([7]) will almost surely go extinct in the subcritical and critical 
case. For the supercritical branching process, the branching process will 
either go extinct or grow exponentially large. In all three cases (subcritical, 
supercritical and critical branching processes) we study a modified process 
which results in an irreducible, aperiodic and positive recurrent Markov 
chain. The stationary distribution of the modified process is the key tool for 
computing E{Ak)- The details of the appropriate modified process differs 
between the three classes of branching process and we therefore consider 
each in turn. We use detailed balance to prove the stationary distribution 
of the modified process which is that same approach that is often used in 
queueing theory to prove insensitivity results, see, for example. Whittle [12j 
and Zachary [13] and references therein. 

2.2 Subcritical branching process: 6 < I 

Theorem 1 1 . 1 1 can be proved in the subcritical case using Zachary |13j . The- 
orem 1. However, it is instructive for studying the supercritical case to 
outline a proof of the result. The branching process will almost surely go 
extinct and we modify ([7]) by regenerating the branching process whenever 
it goes extinct by restarting the branching process with a new individual 
whose lifetime is distributed according to Qi with probability pi. This is 
the regeneration approach introduced in Hernandez- Suarez and Castillo- 
Chavez [7J and extended in Ball and Lyne [2j. That is, for k = ej^„. replace 
k — )• k — ej_„. by k — )• e^ i with probability p/ (1 = 1,2,... ,m). The modified 
process is an irreducible, aperiodic and positive recurrent Markov chain and 
therefore has a unique stationary distribution. Let vTk denote the stationary 
probability of the population being in state k. 



Lemma 2.1 For all k, 



w/iere qij = Pihi.j, C = -5/ log(l - 5) and K = ^™ ^ Y.'jLi Kj- 
For K>1, let Ak = {\^ ■ ET=i E]U hj = K} and 

<Pk= 5Z ^k = CS^'-'/K. (9) 

Proof. To prove the Lemma it is sufficient to show that vTk given by ^ 
satisfies detailed balance. That is, 

y^ TTh/oh.k = TTk y^ Pk,i, (10) 

where ph,i is the rate at which transitions from state h to state 1 take place. 

We start with the right hand side of (jlOp in the case, where k ^ e^^i and 
we do not need to consider the regeneration modifications. 

The possible transitions to state k are from k+ejj— ejj+i (i = 1,2, ... ,m;j 
1,2, ... ,ni — 1), k + ei^n- and k — ej^i. We have the following three equalities. 
For i = 1,2, . . . ,m and j = 1,2, ... ,ni — 1, 

vrk+e,,,-e.,j+i(A:ij + l)7i,i = TTk , ';^ ^ ''^^^ iki,j + lhi,j 
For i = 1,2, . . . ,m. 



ill) 



and 



7rk+e,,„^ (A:i,„, + l)7i,n, = TTkSK ''"-' {ki^n, + l)li,n, 

= ^^^.5Kpi, (12) 



1 A;- 1 

Vrk-e,,i(i^ - l)Pi'5 = VTk-— 7TT^^(-?^ " l)Pi'^ 

(A - Ijd gi,i 

= T^\<ih,iii,i- (13) 



It follows from (fTT HT3]) that 

m rii — l m m 

1=1 j=l i=l i=l 

TTi rii — 1 m m 

= X X vTk/cjj+iTij+i + X TTkSKpi + X vrkA;j,i7j,i 

i=l j=l i=l i=l 

(m rii m 
„■ 1 .■ 1 ,' 1 ^^Lt. 



i=l j=l i=l 



as required. 

For I = 1,2, . . . ,m and k = e^^i, we have that 

m 

i=l 

m m 

= C{6qi^i + P/) X %ri,li,n, = C{6qi^i + 1) ^ Pi 

i=l i=l 

as required. 

Since Eke^,i^!ni^in-=i95VfcM! = (Elti E •=! 9m)'' =J^^^^^- 
lows trivially from [8] and C = — (5/log(l — 6) follows from '^'^=i4'K = 

The above tells us that in equilibrium the modified process spends a 
proportion (pK of its time with K individuals alive. It does not tell us 
directly the mean amount of time, E{Ak), that the branching process has 
K individuals. Each cycle, between regenerations, of the modified process 
corresponds to a single realisation of the branching process. As noted in 
that the mean total size of the subcritical branching process is 1/(1 — 5). 



(14) 



Thus the mean regeneration time, E{T), satisfies 

Y, E{T)^kK = — -, 

K=l 

which gives E{T) = 1/C = - log{l-5)/5 cf. Coroharyini Hence E{Ak) = 
E{T)(f)K = S^~^ / K is the mean time the original subcritical branching 
process with 1 initial ancestor spends with K individuals alive. 

2.3 Supercritical case: 5 > 1 

We modify the approach taken in Section [2.21 to obtain E{Ak) = 'i/{K5) 
when (5 > 1. In this case the expected size of the branching process is infinite 
and there is a non-zero probability of never going extinct. This leads to key 
differences from the subcritical case which need to be resolved. Therefore 
we modify the branching (birth-death) process, defined by dH), to create a 
population process, Vn^ which is restricted to {1,2,..., N} individuals for 
some N &'N. The population process evolves as the branching process ex- 
cept for transitions which lead to or A^ -|- 1 individuals in the branching 
process. In particular, any transition in the branching process which leads 
to A^ -|- 1 individuals (a birth when there are N individuals) leads to the 
population process restarting with 1 individual whose lifetime is distributed 
according to Q. (A catastrophic event killing all A^ individuals in the pop- 
ulation combined with a regeneration event, the birth of a new individual.) 
Also extinction of the branching process (death of the only individual) is re- 
placed by the population process moving to a state with A^ individuals alive. 
(This move does not have a natural interpretation.) The lifetime stages of 
the A'^ individuals will be discussed below. The population process is rather 
different to the regeneration process used for the subcritical case in Section 
12.21 but the resulting Markov chain has a finite state space and is aperiodic 
and irreducible and therefore has a unique stationary distribution. 

Lemma 2.2 For any A^ G N, the population process Vn with stationary 
distribution {vr^} satisfies the following. 

For all k, with 1<K = YT=i TJj=i hj < ^> 

m rii „..^i,j 

^^ = Lr,{K~iy.llll-^, (16) 

• 1 ■ 1 "'i,i' 



where Ln = |l]jii(l/j)| ~ 1/logA^ and w satisfies YT=i YTjLiWij = 

r- m 

li,l+5 = -^ + ^Wl^nill,ni (17) 

*'^ 1=1 

m 

li,j + (5 = -^^^^ ^-^^^— + Y] ^l,ni7l,ni (j = 2, 3, . . . , Hi). (18) 

For 1 < K < iV, /ei Ar = {^^ : YT=i TTjU hj = K] then 

^K= 5^ ^k = Lm/K. (19) 

We comment briefly on the statement of Lemma 12.21 before embarking 
on the proof. Firstly, the form of tt^ given by (J16p is very similar to that 
for vTk given by ([8]). The key difference is that 5 does not feature explicitly 
in (J16p and Wi,j has a more complicated form than qjj = Pi/^ij. For the 
subcritical case gjj is the probability that an alive individual is in the j^^ 
stage of a lifetime distributed according to Qi and Wij plays the same role 
for the supercritical case. In the supercritical case, the branching process 
is growing and we observe that Wi^i > q-i^i and Wi^m < Qi,ni- Fortunately 
it is not necessary to compute w, although for Q ~ Gamma(2, 2) it is 
straightforward to show that wi^i = w and wi^2 = 1 — w, where w = {—6 + 
VS"^ + 8S)/4. Secondly, ([19]) follows trivially from ([16]) and gives us that 
E{Ak) oc 1/K. Therefore after proving Lemma 12.2] we show that the 
constant of proportionality is 1/5. 

Proof of Lemma 12. 2L Fix N £ N. We follow Lemma 12.11 in proving 
the lemma using detailed balance. That is, showing that for all k, vr^ given 
by ([TB]) satisfies 

E ^/'h.k = vr^ E Pk,i. (20) 

We start with the right hand side of (|20p in the case, where k ^ e;^i or 
K = N, that is, we do not need to worry about modifications. 

The possible transitions to state k are from k+ejj— e,j+i (i = 1,2, ... ,m;j 
1,2, . . . ,ni — l), k + ej^„. and k — ej^i. We have the following three equalities 
which are similar to ()llH13p in Lemma l2. II However, there are differences 
due to the different form of tt^ between the subcritical and supercritical case. 



For i = 1,2, . . . ,m and j = 1,2, ... ,ni — 1, 

<e..,-e.,,,, ih, + 1)7.. = ^k 7r^f:^^(fc.. + 1)7. 






For i = 1,2, . . . ,m, 



TT^+e, „. (^i,n, + l)7i,n, = TT^Jf ''"^ (fci,„^ + l)7i,„^ 

m Ui 

= ^k X]^^*'j'^w7w, (22) 



and 



vr^-e.,(i^-l)p.^ = -k(^^(i^-l)P.5 

= ^k • (23) 

The rate of transition out of state k is X^I^i Z]^=i ^i,j^i,j^^^ — J2^i X^?=i ^ijilij^ 
6). Hence, for (j20p to hold, it follows from (J21H23P that we require that 

m rii — l I m rii \ m ^^ U X 

E E ^^.^>^|;J]^ + E E ^.. E --^7... + E ^ 
1=1 .=1 '■' \ j=i j=i / 1=1 1=1 ' 

■m ?ii 

= EE^m(7m+'5)- (24) 

i=i j=i 

Equating the coefficients of the kij terms we get w solving ()17p and (jlSp . 
Therefore we need to check that this choice of w is also consistent with the 
boundary cases. 

For k such that K = N, we have that 

m m m rii -^j 

E<.,^'.«^^^/(^^/^)} = Y.^NWl,,,^l,n,LN{N-l)\\[\{f^X — 
1=1 1=1 1=1 j = l '-' 

m 



X](^^ + l)7i,n,^k+e,,„,, (25) 



1=1 
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where, with an abuse of notation, we take vr^^ to satisfy ()16p with K = 
N + 1. That is, the transitions from a single individual to N individuals 
in the modified process mimic the transitions in the branching process from 
A^ + 1 individuals to N individuals. All other transitions in and out of state 
k are identical to the branching process and it is straightforward to verify 
that (I20|) holds. 

For k = ei 1, we have that the transitions into state ei i are from e; i + 
Gi^m and from any k such that K = N. Using (j)^ = X^^g^ ^k^ ~ L^/N, 
the rate of entry into state e;_i is 

PiAf<5(/)jv + 2^7i,n.vr^,,+e,,„^ = PiLnS + l^-/i,n, T • (26) 

1=1 ' ' i=l 

The rate of exit from state e^ i is (7/,i + (^). Hence, from (|26|) we require that 

m 
PiLn6 + ^ Ji,n,L]yWl^iWi^n, = {ll,l + 5)LnWi^i. (27) 

i=l 

However, (f27|) is equivalent to (fT7|) . (Simply multiply both sides of (fT7|) by 
LnWi,!-) Therefore completing the proof that vr^ is indeed the stationary 
distribution of the modified population process. D 

We now show that E[Ax) = 1/{5K) by studying the total amount of 
time, T/v, that a supercritical branching process spends with between 1 and 
A^ individuals alive. In particular we study the asymptotic behaviour of 
Tat as A^ — )• oo. The supercritical branching process either goes extinct with 
probability, z, or explodes {X{u) — )• oo as u — )• oo) with probability, 1 — z. It 
is helpful to consider these two cases separately with Ext denoting the event 
that the branching process goes extinction. Firstly, it follows from Cham- 
pagnat et al. [3j, Proposition 2.1, that z = 1 — r]/6 , (see also Lambert [9], 
Section 5), where r/ is the Malthusian parameter of the branching process. 
Moreover, Lambert [9], Proposition 5.7 gives the law of the supercritical 
process conditional upon extinction which is a subcritical branching process 
with a modified lifetime distribution and birth rate 5 — tj. Thus it is trivial 
to show that E{Ti\f\F,xt) = 0{1). Conditional upon non-extinction, Cham- 
pagnat et al. [3j, Proposition 2.2 states, exp{—r]u)X{u) ^-> y as u — )• oo, 
where Y is exponential random variable with mean dependent upon r/ and Q. 
(This is proved in Lambert [9] and is a special case of Nerman jllj . Theorem 
5.4.) Let T^ = in.m{u; X (u) = A^ + 1} and T^ = m.ax{u;X{u) = N} denote 
the first time X{u) leaves {1, 2, . . . , N} and the last time X{u) belongs to 
{1, 2, . . . , A^} with T^ < Tj\f < T^. Then conditional upon non-extinction, 

11 



it is straightforward to show that T^ x 77/ log N, T^ x rj/ log A^|Ext —^ 1 as 
A^ — ;■ 00 and consequently that T/v x t]/ log N{= TAr)|Ext —A- 1 as N ^ 00. 
Given the structure of the branching process it is straightforward to show 
that there exists < g < 1, such that for all m G N, P{Tn > m) < q^, and 
hence, that -E(TXr|Ext ) — )• 1 as A^ — )• 00. Therefore 

E{Tm) = P(Ext)^(T7v|Ext) + P(Ext^)-E(rjv|Ext^) 

floe TV 
= zO(l) + (l-z)<^^^^^^ + o(logiV) 

■n logN ,, ^^, 
= -jx ^^ + o{logN) 
rj 

= -logiV + o(logiV). (28) 

Since the above holds for all A^, we have that 

EiAK) = lim cP^E{Tn) = hm ^^JlogA^ = ^ (29) 

Af-s-oo N^oo K oK 

as required. 

2.4 Critical case: 6 = 1 

For the critical case the expected duration of the branching process is infinite 
and the probability of non-extinction is 0, which makes studying this case 
particularly difficult. The modified population process defined in Section 
12.31 can still be constructed in this case and Lemma 12.21 still holds with 
1i,j{— Villi,]) = ^i,r Thus it is straightforward to show that E{Ak) oc 1/K 
with lim^ii 6 /K = lim^^i 1/{6K) = 1/K. However, it is difficult to 
adapt the approach taken in Section [2.31 to show that E(Ak) = 1/K but to 
prove the result it suffices to fix Q and to show that E{Ai) = 1. 

We study the time that the branching process spends with one individ- 
ual alive, which we term the local time process (of the branching process 
with one individual). Suppose that Q satisfies (0). We say that the local 
time process is in state (i, j) if the individual is in the j*'' stage of lifetime 
Qi. The local time process regenerates with a new individual in state (i, 1) 
with probability pi {i = 1,2,..., m) if the branching process goes extinct. 
Also since the branching process is critical it will almost surely go extinct. 
Therefore the branching process started from one individual will eventually 
return to one individual, although the mean waiting time is infinite. The 
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transitions from state {i,j) in the local time process are given by:- 

{(a, 6) rate r(ij)_(„_fe) 

{i,j + 1) rate jij if (j = 1, 2, . . . , n^ - 1) (30) 

(a, 1) Tate -fijpa if j = Ui 

where X^(a m '^(jj),(a,b) = 1 ^^d ?'(ij),(a,6) is the probability that following a 
birth to an individual in state {i,j) the next time the branching process 
returns to one individual, the sole individual will be in state (a, 6). The 
matrix R = {rtij\fa,b)) is difficult to compute but we know the stationary 
distribution for the local time process from the stationary distribution of 
the branching process. Hence the probability of being in state {i,j) is qij = 
Pi/lij- Therefore following Ball and Lyne [2], Section 2.2.3, we note that 
the overall intensity of regeneration of the local time process (extinction of 
the branching process) is Y2^i Qi,ni'yi,ni = Y^^iPi = 1- Thus the mean time 
between regenerations is E{Ai) = 1~^ = 1 as required. 
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